Transient localization in crystalline organic semiconductors. 
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A relation derived from the Kubo formula shows that optical conductivity measurements below the 
gap frequency in doped semiconductors can be used to probe directly the time-dependent quantum 
dynamics of charge carriers. This allows to extract fundamental quantities such as the elastic and 
inelastic scattering rates, as well as the localization length in disordered systems. When applied 
to crystalline organic semiconductors, an incipient electron localization caused by large dynamical 
lattice disorder is unveiled, implying a breakdown of semiclassical transport. 

PACS numbers: 



Introduction. "Bad" conductors are systems present- 
ing a breakdown of the semiclassical Bloch-Boltzmann 
description of electronic transport. Known examples are 
found in various classes of materials such as disordered 
systems [l| , transition metal compounds Q , alkali-doped 
and charge-transfer organic metals [1] and quasi- crystals 
0, HI- In all these cases, the electron motion is so much 
slowed down (by disorder, electronic correlations, pola- 
ronic effects, or by structural constraints) that the semi- 
classical assumption of well-defined wave-packets under- 
going rare scattering events is not valid. 

An analogous situation is encountered in crystalline 
organic semiconductors. There, a fundamental unsettled 
question is whether the mechanism of charge transport 
can be ultimately understood from the point of view of 
band electrons alone, as suggested by the "band-like" 
temperature dependence of the measured mobility. In 
these materials, the inherently large thermal molecu- 
lar motions act as strong electron scatterers, leading to 
apparent electron mean-free-paths comparable or even 
smaller than the intermolecular distances. While a gen- 
eralization of semiclassical transport theory has been re- 
cently proposed in Ref. @ to deal with this situation, it 
might well be that the very nature of semiclassical trans- 
port is insufficient to appropriately describe the charge 
transport mechanism in these materials. Indeed, recent 
numerical studies have suggested a radically different 
point of view, based on a form of electron localization 
due to the dynamical disorder caused by the thermal 
molecular motion [7], |8J. Accordingly, a theory of elec- 
tron transport in organic semiconductors would require 
a proper account of quantum corrections to the electron 
dynamics, not included in semiclassical treatments. 

In this paper we tackle this problem by expressing the 
Kubo formula as a relation between the optical conduc- 
tivity and the time-resolved quantum dynamics of elec- 
trons. We first consider a microscopic model with dy- 
namical lattice disorder that illustrates the character- 
istic behavior of the quantum diffusion. A relaxation 



time approximation is then introduced that treats the ef- 
fect of inelastic scattering by low-frequency lattice vibra- 
tions in an intuitive way. Finally these concepts are used 
in the interpretation of experimental data in crystalline 
rubrene, providing evidence for localization effects. 

Formalism. The quantum diffusion of electrons in 
a given spatial direction can be measured via their 
quantum-mechanical spread 



AX 2 (t) = ([X(t)-X(0)} 2 ), 



(1) 



where X(t) = Yli=i%i(t) 1S the total position opera- 
tor of N electrons in the Heisenberg representation and 
(••■) = Trfe - ' 9 ^ (. . .)]/Z denotes the thermodynamic av- 
erage. AX 2 (t) is directly related to the symmetrized 
self-correlation function of the velocity operator Vx (t) — 

C(t) = (Vx(t)Vx(0) + Vx(0)Vx(t)), via 0| 



dAX 2 {t) 
dt 



C{t')dt'. 



(2) 



On the other hand, the Kubo formula expresses the dis- 
sipative part of the optical conductivity as 



o-(cj) 



— — Re 

vTioo ,/q 



^c(0)])di, (3) 



where v is the volume of the system. This can be ex- 
actly related to the symmetrized C(t) of Eq. ©, and 
therefore, to the quantum diffusion AX 2 (t). Replacing 
the commutator in Eq. ([3]) with the anticommutator of 
Eq. @ can be absorbed into a detailed-balance prefac- 
tor, yielding: 



cr(w) = - 



e 2 w 2 tanh(/??kj/2) 
v Tiuj 



Re 



l AX 2 (t)dt (4) 



with (3 = 1/ksT (see Refs|9|, llfj for related aspects). 
The above Eq. ([4| is a restatement of the Kubo for- 
mula, identifying the time-dependent quantum diffusion 
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FIG. 1: a) Quantum spread obtained from the microscopic 
model Eq. §6§ for static (grey, dashed) as well as dynami- 
cal disorder (red, full lines: from bottom to top, huo/J — 
0.01,0.0435,0.1). Times are in units of h/J. b) The corre- 
sponding instantaneous diffusivity D(t) = (dAx 2 /dt)/2. The 
inset shows the optical conductivity obtained via Eq. @. 



as the physical quantity that is dual to the optical ab- 
sorption in the frequency domain. For independent non- 
degenerate electrons, the formalism presented above ac- 
quires an intuitive meaning in terms of the quantum 
spread of the electronic wavefunctions, as in this case 
one has AX 2 (t) = NAx 2 (t), with Ax 2 (t) referring to 
each individual particle. Eq. (0} can be inverted to give 



Ax 2 (t) = -^Re [ 
ne 2 Jq 



a(ui)/n 



u> tanh(/3?iu;/2 



■duj. (5) 



with n — N/v the electron density. 

Microscopic description of lattice disorder. To set the 
concepts clear, we apply the quantum diffusion formalism 
to the following Hamiltonian [(J 0] , 



H 



t+l)] {ct C i+l+ C t+l c i)+ H ph- (6) 



In this model, electrons in a one-dimensional tight- 
binding band have their inter-molecular transfer inte- 
grals modulated by molecular vibrations of frequency 
ujq, represented by H ph = M "o u i + p±__ g q flg 
captures the effects of strong dynamical lattice disorder 
that seem to be crucial in crystalline organic semicon- 
ductors. The phonon frequency wo is small due to the 
large molecular weight, so that the lattice fluctuations 
can be treated classically (fc^T > Hloq). Their coupling 



to electrons is governed by the dimensionless parameter 
A = a 2 J/{2Muj 2 ). 

To calculate the electron diffusion Ax 2 (t) in the pres- 
ence of lattice dynamics we employ mixed quantum- 
classical simulations based on the Ehrenfest coupled 
equations on a 1024-site chain. We solve the 
Schrodinger equation and average over up to 12800 initial 
conditions, with the initial displacements, Ui, obeying a 
thermal distribution P(w 4 ) oc exp(-MwQU^/2fc_B? 1 ). The 
case of a frozen disordered lattice is treated by averaging 
over the same set of disorder realizations (see Ref. [6]). 
In the latter case, the results are cross-checked via an ex- 
act diagonalization of the electronic problem on 256 sites, 
allowing for a direct verification of Eqs. (|4|) and ([5]). 

Fig. Q] shows the electron spread Ax(t) = Ax' 



dAx-'(t) 
2dt 



'(*) 

for 



and the time-dependent diffusivity D(t) = 
a representative choice of microscopic parameters: A 
0.25, J = O.lleV, T = 0.235J = 300K, and different 
values of ujq. The dashed line in Fig. QJ, is the result 
for static disorder, ujq = 0, showing a finite localization 
length L = 3.6a. The corresponding D{t) in Fig. [T]d 
increases at short times in the ballistic regime, then ex- 
hibits oscillations. At subsequent times the oscillations 
are damped and D(t) steadily decreases and vanishes. 
Following Eq. (2), 2(dD/dt) = C(t) is precisely the ve- 
locity correlation function: a negative slope is therefore 
signalling the occurrence backscattering underlying the 
phenomenon of Anderson localization (in a classical pic- 
ture, C(t) < implies that the velocity at time t is oppo- 
site to its value at time t — 0). This occurs, as expected, 
at times greater than the elastic scattering time [l[ , which 
is given by r e ; = (nXT)^ 1 = 5.4 in the present units 

The results in the presence of lattice dynamics (ujq ^ 0) 
closely follow the localized behavior at short and inter- 
mediate times. However, upon reaching the timescale of 
lattice vibrations, 1/wo (indicated by arrows), localiza- 
tion is destroyed and Ax(t) starts increasing indefinitely. 
The existence of a transient localization phenomenon at 
times Tei ^ t ^ l/wo is one of the main results of this 
work. It indicates that the electronic transport mech- 
anism is markedly non-semiclassical, the final outcome 
being determined by the characteristic timescale of lat- 
tice disorder. 

As a side remark, Fig. [TJ) illustrates a fundamental 
drawback of the Ehrenfest method, that makes it inad- 
equate to determine asymptotically the electron diffu- 
sion: the diffusivity does not apparently tend to a con- 
stant value but rather exhibits an upward drift at long 
times. The total energy of the system (not shown) is 
conserved in the simulation with a relative precision of 
2 • 10~ 7 , which rules out possible integration errors. The 
origin of this spurious phenomenon rather lies in the fact 
that the Ehrenfest equations do not properly conserve the 
Maxwell-Boltzmann statistical distribution 11] : the re- 
peated action of the lattice vibrations (an external, time- 
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dependent potential) eventually drives the electrons to an 
arbitrarily high effective temperature, a fact that could 
be at the origin of the T~ 2 temperature dependence of 
the mobility obtained by this method 0, fl2l. Il3| . 

The inset of Fig. [Tp shows the optical conductivity, 
cr(ui), calculated by applying Eq. to the data of 

FigQJ,, neglecting the spurious superdiffusive behavior 
at long times. The result of the static disorder problem 
(luq = 0) is shown for reference (grey, dashed). We see 
that the dynamical nature of the lattice only modifies the 
low frequency region of the spectrum, u < loq. It does not 
affect substantially the localization peak at uji oc ~ 0.4J 
as long as ujq <C w; c, nor the absorption band at higher 
frequencies. 

Relaxation time approximation. To understand how 
localization features can actually coexist with a diffusive 
behavior at long times, we now implement the relaxation 
time approximation (RTA) as a simple scheme bridging 
between localization and diffusion. The idea underlying 
the RTA is to express the dynamical properties of the ac- 
tual system in terms of those of a suitably defined refer- 
ence system, from which it decays over time. Specifically, 
defining Cq as a reference velocity correlation function, 
the relation 



C RTA {t) = C {t)e- t ^ 



(7) 



describes the damping of velocity correlations caused by 
relaxation processes with a characteristic time r 0, HI, 



14j . In the semiclassical theory of electron transport, one 
starts from a perfectly periodic crystal and describes via 
Eq. ([7]) the momentum relaxation due to the scattering 
of Bloch states. In that case Co = 2i> 



avy 



is a constant 



(twice the thermal average of the squared band velocity) 



-t/r 



and the resulting diffusivity, Drta^) 
is a monotonically increasing function of time. 

One can alternatively take a localized system with 
static disorder as the reference state, as suggested by 
Fig. [TJ At times shorter than the typical timescale of 
the lattice motion, Tj n ~ 1/woj the molecular lattice ap- 
pears to the moving electrons as an essentially frozen, 
disordered landscape. In this case, nothing prevents the 
buildup of quantum interferences that lie at the origin 
of Anderson localization. The velocity correlation func- 
tion C(t) then initially coincides with that of a system 
with static disorder, Co(t). However, quantum interfer- 
ences are destroyed at longer times because, due to the 
lattice dynamics, the electrons encounter different disor- 
der landscapes when moving in the forward and back- 
ward directions [lj]: Eq. ([7]) is the simplest form that 
is able to capture such decay process. From Eqs. ([2]) 
and (|4]) it is easy to see that, starting from a local- 
ized system (i.e. one with a vanishing diffusion constant, 
/ °° Co(t)dt — 0), Eq. (UJ restores a finite diffusion con- 
stant, Drta = Lg(Ti n ) I \2rin) which is analo gou s to the 
Thoulcss diffusivity of Anderson insulators 15]. This 




FIG. 2: Color online. Time dependent electron diffusivity 
D(t) extracted via Eq. ((4} from the experimental optical 
conductivity of Ref. [Hj in the direction of highest conduction. 
The absolute value is fixed by the measured mobility \i ~ 
7cm 2 /Vs. The inset shows the same quantity as a function 



of the instantaneous electron spread, 
weak localization extrapolation. 



The dashed line is the 



value is essentially equal to the diffusivity of the refer- 
ence localized system at a time t w Tj n . The quantity 
^o( T m) = / e~*/ Ti " AxQ(t)dt/r in represents the typical 
electron spread achieved at a time Tj n , before diffusion 
sets back in. It therefore acquires the meaning of a tran- 
sient localization length. The emerging physical picture 
is that of electrons prone to localization, but that can 
take advantage of the lattice motion to diffuse freely over 
a distance Lo(Ti n ), with a trial rate 1/Ti n . 

From Eq. ((U and Eq. ((Jj we obtain a mobility 



e Ll{T in ) 



(8) 



Although a systematic study of the electron mobility of 
organic semiconductors is beyond the scope of this work, 
we note here that under quite general assumptions, Eq. 
© implies a power-law temperature dependence, /i ~ 
T~", even though the microscopic transport mechanism 
is far from conventional band transport (the exponent a 
depends on how the transient localization length varies 
with the thermal lattice disorder). 

Real time dynamics from experiment. We now show 
how optical conductivity experiments can provide direct 
information on the relevant time and length scales of the 
problem. Fig. [2] reports the instantaneous diffusivity, 
D(t), in the direction of highest conduction of rubrene, 
obtained via Eq. f5| by direct integration of the data of 
Ref. [IH (analogous results are obtained from Ref. Klh . 
D(t) increases first, reaches a maximum and then de- 
creases by a factor of 3 before stabilizing to a constant 
value. The shape of the diffusivity curve is remarkably 
similar to the theoretical result of Fig UJ and hardly com- 
patible with the semiclassical picture discussed after Eq. 
(UJ in which D(t) increases and directly saturates. This 



4 



suggests that, as in the model calculation, the decrease 
of D{t) is due to localization effects occurring at times 
shorter than the lattice dynamics. An elastic scatter- 
ing time of the order of r e i « 10 _14 s can be tentatively 
identified with the region of negative slope in Fig. [2] Ac- 
cording to the arguments given above, a diffusive regime 
\D(t) =constant] sets up at timescales beyond the inelas- 
tic scattering time. This is what Fig. [2] reveals, with 

T in « 5 • 1CT 14 S. 

The inset shows a plot of D(t) as a function of Ax(t) 
and gives access to the relevant length scales. We find 
for the elastic and inelastic mean-free paths £ e i « a and 
£ in « 3a respectively, with a = 7.2A the intermolecu- 
lar distance. In addition this inset suggests that with- 
out inelastic scattering, i.e. for a fixed set of disordered 
molecular positions, the diffusivity D(t) would extrapo- 
late to zero at a localization length of the order of 3 — 5 
intermolecular distances. Note that the linear extrapo- 
lation of diffusivity with the Log of the length is a stan- 
dard approximation for two dimensional systems [l[ . Yet 
rubrene should be considered as intermediate between 
one and two due to its highly anisotropic character [lij . 
For a one-dimensional system localization sets in more 
efficiently, therefore the above extrapolation should be 
an upper bound to the true localization length. 

We finally show in Fig. [3] how the RTA can be used 
to extract quantitative microscopic information from the 
optical data. One first constructs an ansatz for the refer- 
ence conductivity <To( w ) representing the ideal case with 
frozen disorder, i.e. no inelastic scattering. This can 
be done starting from the experimental optical absorp- 
tion (inset: red, full line), by enforcing the condition 
cjq{uj — > 0) = appropriate for a localized system (grey, 
dashed line). The RTA result (black, dash-dotted) is 
then obtained by applying Eq. ([7} to fit the experi- 
mental curve. The optical conductivity and the result- 
ing diffusivity both nicely agree with the experimental 
data. The fitting procedure yields T in = 5.1 • 10~ 14 s, 
corresponding to a frequency u>i n — 104cm , consis- 
tent with the relevant intermolecular phonon frequencies 
in rubrene [l|| [2(|. From the same fit, the estimated 
transient localization length in the direction of highest 
conduction is Lo(7m) ~ 2a [the localization length of the 
static ansatz is L$(t — > oo) ~ 3a]. The present analysis 
shows that the finite frequency absorption peak observed 
in rubrene ill 17| should be ascribed to transient local- 
ization effects — i.e. occurring before the dynamics of the 
lattice set in — and constitute a signature of an uncon- 
ventional transport mechanism. 

Concluding remarks. The relation between the quan- 
tum dynamics of electrons and the optical conductivity 
that stems from the Kubo formula, appears to be a pow- 
erful tool to analyze the charge dynamics in semicon- 
ductors with unconventional transport properties. When 
applied to experimental data on crystalline organic semi- 
conductors, it provides evidence for the role played by 
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FIG. 3: Color online. The experimental diffusivity of Fig. [2] 
(red, continuous) is compared with the result of the relaxation 
time approximation (RTA) (black, dash-dotted) and the local- 
ization ansatz defined in the text (grey, dashed). The inset 
shows the experimental optical conductivity of Ref. 03 to- 
gether with the RTA and the localization ansatz (arb. units). 



localization phenomena in the charge transport mecha- 
nism. The scenario emerging from the above analysis 
is indicative of a prominent role of the dynamical lattice 
disorder, which is supported by a microscopic calculation 
on a one-dimensional model. 
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